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S»J>tH\RY 


Th« crack growth analysis program davsloped for the George 
C. Marshall Space Flight Center by Del West Associates, Inc. 
is described. The technical approaches used within the 
computer program are presented and the input and output 
formats and options are described. Example data and example 
problems are included. 
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INTROOUCTIOM 


In order to include the important consideration of structural 
failure du't to the presence of flaws and crack-like defects in 
aerospace hardware; it is necessary to have a computer program 
capable of performing crack growth analysis that is easy to use 
and generally applicable. The need for a computer program (as 
opposed to simple hand calculations) arises from the complexity of 
growth descriptions required for crack growth analysis of real 
materials in complex structure under a variety of loading and en- 
vironmental conditions. The MSFC crack growth computer program 
developed by Del West was designed to meet this need. 

The MSFC crack growth computer program calculates crack growth 
for part through cracks, through the thickness cracks and cracks 
which are tr^.nsitioning from part through cracks to through the 
thickness cracks. The computer program has been written to be 
flexible in its operation and to be easily adapted and changed 
as fracture mechanics technology changes and/or the design usage 
of the program changes. 

The computer program Js essentially an integration routine which 
calculates crack growth from an in initial defect size and ter- 
minates calculation when the crack is sufficiently large for a 
critical condition (instability or rapid growth) to be reached. 

In addition, if a design life is not met for a particular structure, 
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the program has the capability of varying the thickness of the 
structure so as to establish the thickness which will meet the 
design requirements. 

During the period when a crack Is a part through crack, crack 
growth in iLe depth and surface directions may be different due to 
variation* in stress intensity factors and/or directional depend- 
ence of material properties. The MSPC computer program considers 
both of these effects and hence incorporates realistic crack shape 
changes. During the period when a crack is transitioning from a 
part through crack to a through the thickness crack, the crack 
lengths on the backside and the frontside are different. The MSPC 
computer program tracks the growth of these two dimensions separately; 
evaluating the stress intensity factors at each surface until these 
din.*r;.sions are the same and the crack has completed its transition 
to a through the thickness crack. 

The computer program allows two different methods of load input. 

For each step in the loading block, the user specifies either; 

(1 ) Maximum Stress, Minimum Stress, N/inber of Cycles or (2) Maximum 
Stress, Stress Ratio, Number of Cycles. It should be noted that if 
crack growth mechanisms other than fatigue are being considered 
(e.g., static stress corrosion) the appropriate rate variable can 
be used instead of cycles (e.g., time at load) in conjunction with 
appropriate material constants as described below to perform a 
wide range of phenomenological studies. 
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Th« us* of s limit losd (a load which may b* higher than any load in 
th* actual spectrum) to determine the end of design life is a 
common practice. The MSFC computer program has therefor* been 
written to consider a separate limit load (apart from those in 
the spectrum) and to determine when it causes failure, ttowever, 
after failure due to limit load occurs, the crack growth cal- 
culation continues. The limit load failure information is included 
in th* output. 

The crack growth rate material properties may presently be input 
into tho program in any of three formats! (1) Paris equation with 
upper and lower cutoffs in stress intensity factor; (2) Forman 
equation with upper and lower cutoffs in stress intensity factor J 
(3) Collipriest-Ehret equation with additional upper and lower cut- 
offs in s*;Yess intensity factor. An important feature of the 
material property description is that different materials properties 
(crack growth equations, fracture properties, yield stress, etc.) 
may be designated for each step in the loading spectrum. Thus varying 
temperatures and environments may be considered. 

it 

The MSPC crack growth computer program has tho capability of 
utilising any one of three crack growth retardation models. Of 
course, the effects of retardation on crack growth will not be 
considered if the user does not request it. The three models 
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prasently available are: (1) Willanborg; (2) Wheeler} (3) 

Grumman Cloaure Model . 

The module which performs stress intensity calculations currently 
includes stress intensity equations fort (1) Part through cracks 
with front and back surface correction factors; (3) Transition 
cracks; (3) Center cracked panels with a finite width correction 
factor; (4) Compact test specimens. 

As many runs (each with varying input conditions) as desired may 
be stacked. As additional runs are made, only that section of 
data which is changed (i.e., loads, material* properties, or geometry) 
need be reentered. Output for each run nay be controlled by the 
user. Minimum output fox each run consists of information on input 
data and failure (crack lengths, cycles, etc.) as well as crack 
lengths, stress intensity factors, and crack growth rates for the 
first and last cycle of each stress level in the first load block 
applied as a part through crack, transitional crack or through 
crack. Additional information (crack lengths, stress intensity 

* A 

factors, and crack growth races) for particular blocks and loading 
steps may be requested by the user. 

A flow chart showing all subroutines is presented in Figure 1. 
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Subroutln* Function 

M/VIN Reads Input, Sequences Runs, Performs Iteration on 

Thickness, Calls Appropriate Crack Growth Module. 

PTC GRW . Calculates Crack Growth for a Part Through Crack 

TRANS Calculates Crack Growth for a Transitional Crack 

TC CROW Calculates Crack Growth for a Through the Thickness Crack 

KANAL Evaluates all Stress Intensity Factors. 

DAMAGE Calculates Crack Growth Rates 

RETARD* Modifies Input to DAMAGE to Account for Retardation 

Effects. 

Figure 1. Overall Flow Diagram 

MSFC Crack Growth Analysis Program 
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TtCHNICAL APPWJAai 


Th« ••••nc* of tha crack growth analysis procodurss consists of t 

1) Considering each loading step in a load block in turn. 

2) Evaluating stress intensity factors, using the stresses from the 
step under consideration. 

3) Using these stress intensity factors (and previous loading history 
if retardation is considered) to calculate crack growth rata. 

4) Consider a small amount of growth (~15 of current crack sise) 
and calculate the number of cycles it takes to grow that amount. 

If that amount exceeds the lumber of cycles not yet consumed in the 
step then only those remaining cycles are used and a corresponding 
crack growth increment is calculated. 

5) Crack lengths are incremented, cycle count is incremented. 

6) This process is continued until all cycles in the step are con- 
sidered. The next step is then called. At the end of a block the 
first step is called again. 

7) The calculation ends when: 

i\) The critical stress intensity (s.‘ther at the surface or at 

the depth of a crack) is exce ded. 

•■8 

b) There is no crack growth ( <10* in.) for an entire block. 

c) The crack growth rate goes to infinity (when using the Forman 
equation for crack growth rate). 

d) The maximum number of blocks is exceeded. 

8) All input and all output data are in units conq^atible with Kips 
and inches, (e.g., Ksi, Kai ^in. and in/cycle.) 
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Th« aubroutincs PTCX2RW, TRANS, and TCGROW calculata tha cr%ck growth 
incrananta, raturn to main for inforaatlon on tha naxt loading atap, 
conaidar whan to and tha calculation and tranafar to aach othar 
(PTCQRW -* TRANS TCCROW) aa raquirad. For a part through crack 
PTCORW parforma thaaa functiona until TRANS ia callad. TRANS ia called 
whan tha crack depth aquala tha plate thickneaa. TRANS parforma thaaa 
functiona while tha crack ia tranaitioning to a through crack and calla 
TCGROW whan tha back aurfaca length excaeda 9591 of the front aurfaca 
length. TCGROW parforma these functions whan tha crack ic a through 
crack and may be callad by TRANS or in those cases whan a through crack 
is considered initially it is called from MAIN initially. 

Tha aquations in KANAL, DANVtGlv AND RETARD are, of course, tha heart of 
tha crack growth analysis program, These are described below. 

KANAL 


KANAL is a subroutine which returns factors, which when multipied by 
the appropriate loading term yields stress intensity factors. Thus the 
loading input must be conpatable with the crack configuration con- 
sidered. For the configuration currently in the program the corres- 
ponding name (iCTYPO) and required load description are aa follows: 

KTVro Configuration Load Description 

I Center Cracked Gross Stress 

Tension Panel (^i) 
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ICTYPO 

Configuration 

Load Description 

2 

ASTM B .399 

Pin Force 


Compact Specimen 

(Kips) 

3 

Part Through 

Cross Stress 


Crack 

(Ksi) 

4 

Transition Crack 

Gross Stress 



(Ksi) 

The equation used for 

the center cracked panel is 

taken from Feddorson 


(Ref. 1 ) and is t 

K ■ Sec ^ 

where the diraonsions are given in Figure 2. 

The equatir>n used for the compact specimen is taken from Ref. 2 



For the part through crack the stress intensity is evaluated at the 
surface and at the depth. In both computations a term "Q" is used. 

"0" is composed of an elliptical integral, ^(a/2C) which accounts for 
some effects of crack shape and Q also contains a plasticity correction 
factor. An approximation developed by Rawe (Ref. 3 ) for the pertinent 
elliptic integral is used in the computer program. 

^ ^ ■ 1. +4.593 (a/2C)^*^* 


9 



Compact ^>ecimen Part Through Crack 


Figure 2 • Crack Configuration 





The plasticity correction factor is that proposed by Irwin(Ref. 4 ) 
except that the range in stress is used since the stress range and 
the reversed plastic sone it produces is significant in fatigue. 

The equation for Q is 
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Where .Ir/ is the gross stress range ando^^ is the ,2% offset yield 
stress. For the stress intensity at the front surface the equation 
used is 


K (Surface) « (1.12 ♦ .11 a/C)»;V^ § 


Where the expression in parenthesis is a "front surface correction 
factor." This factor is simply a linear variarion with shape between 
the two available solutions (e.g., 1.12 for a/c » 0 and 1.23 for 
a/2 ■ 1, Ref . 5 and 6 respectively.) The expression for the stress 
intensity at the depth is taken from the work of Shah and Kobayashi 
(Ref. 7 ) with K (depth) ■ n ^front 

Their expression for * 

Mfront » 1 ♦ .2 (I ” a/2C)^ 

s 

Their backface correction factor has been put in an equation form by 
Collipriest and Ehret (Ref. 8 ) and 1st 


^ack 


= 1 + 


j^[o.089 (f).0.2315(i)' 

- 0. 3873 (f)’ * 5. 28 -9.11 (f)' * 5. 233 (i)‘j 

X [l. 109 - 9. 142 (-^) 4 41. 56 - 86. 55 (-^)’ 
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Th« stress intensity equations used for a crack during transition 
are approximate expressions developed by R. M. Ehret which predict 
observed crack growth behavior. The stress intensity for the front 
surface if taken to bet (See Figure 2 ) 


K. ^ m ^nC Sec 
front ^ 


2W 


and the stress intensity for the back surface is taken to be 
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rwMAGe 


The subroutine DAMAGE currently contains three equations for cal- 
culating crack growth rate. In all cases the independent variables 
are the effective stress intensity factor, KE, and the effc«.tive 
stress ratio, RE. IVhen retardation is not used KE is simply the 
stress intensity range (KM\X - KMIN^ and RE is siiiq>ly the stress 
ratio. When retardation is used, KE and RE are calculated in 
RETARD. 


The following equations all contain material property constants 
designated by D(NC,1,J). indicates whether the surface (NC *1) 

or depth (NC ■ 2) is being considered. 1 i« the constant number in 
that equation and indicates the order of the constant on input type 
10 cards. J Is the material type number. In order to call out the 
proper equation the corresponding equation name (NEQ) must be called 
out on card typo 9 columns 11 - 14. Note that each material type 
could use a different NEQ. 

NEQ EQUATION 

1 Collipriest-Ehret 

2 Paris 

3 Forman 
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The CoIllpriost'Ehret equation Is 


dn 


Cl EXP 


C2 


TANK 


-I 


/ ln(KE^/n-RE) D (NC.3.J) D (NC./«..n \ [ 

I In (l-RE) D (NC,A,J)/D (NC.3.J) J [ 

C,l,J) ^(NC,3,J) D(NC,«,J)j 

\ D(NC,4,JV 


D(NC^2>J) 


C 2 ■ In 


D(NC,2,J) 

2 


where 


D(NC,l,J) 

0(NC,2,J) 

D(NC,3,J) 

D(NC,4.J) 


Crack growth rate coefficient 

DlmcnslonleDU Coefficient relating to 
midrange slope 

Critical stress Intensity (upper asymptote) 

Threshold stress Intensity range (lower 
asymptote) 


The Paris equation is 


da 

dn 


D(NC,l,J) 


ypD(NC,2,J) 


The Forman equation Is 


^ ^ D(NC,2,J) 

da _ D(NC.l.J) Ki: 

dn “ (1-RE) D(NC,3,J)-KE 


where 

D(NC,1,J) 

D(NC,2,J) 

D(NC,3,J) 


Crack growth rate coefficient 

Crack growth rate exponent 

Critical stress Intensity (upper asymptote) 



RETARD 


The subroutine RETARD currently contains three retardation models. 

In each of these a crack tip plastic sone r^ is calculated accord- 
ing to the equation 

2 

I V K max 

% -(2T)?7- 

where IPs is a constant depending on the degree of plane stress 
versus plane strain For plane stress Ps ■ 1. For plane strain 
Ps ■ 3. 

The following is the retardation equation number for each retardation 
model. (Input for Columns IS - 18, Card Type 9) 

NRET 

1 Willcnborg Model (Ref. 9 ) 

2 Wheeler Model (Ref. lo ) 

3 Grumman Closure Model (Ref. 11 ) 

The Willenborg retardation model calls for no constants other than 
Pz. The only material property input for this model is therefore 

CR(NC,1,J) » Pz 

The "Wheeler" model in this computer program is actually a variation 
of the model originally presented by Wheeler. Wheeler used a modifi- 
cation to the crack growth rate to produce a retardation effect and 
we have used a modification to the dependent variable KE. If the 


15 


Pari* *quation is u**d for crack growth rat* th* "Wh**l*r" mod*l 
in this computar program i* idantical to th* modal praaantad in 
Raf, 10. Th* input matarial propartia* for this modal ara 

CR(NC,l,J) ■ Pa 
CR(NC,2,J) ■ m/n 

whar* m is iJantical to the "m" used in Ref. 10 and n is th* 
axponont in the Paris equation (i,e., if the "Wheeler n" were 
S and the Paris "n" 4, the input value for CR(NC,2,J) is 1.25) 

Th* details of the Grumman Closure model are too complex to be 
described here. The input is described below 


Input Quantity 

Name in Re 

CR(NC,1.J) 

Pa 

CR(NC,2,J) 


CR(NC,3,J) 

O 

CR(NC,4,J) 

P 

CR(NC,5,J) 

NSAT 

CR(NC,6,J) 


CR(NC,7,J) 

B 
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INPUT 


Th«r« tktm 13 distinct data input card formats. Thsss ars dsscribsd 
bslowt 


CARD FIELD 
TYPE 

1 1-40 

2 1-4 

2 5-10 


2 11-14 


2 15-18 


3 1-10 


3 10-14 


NOMENCLATURE 


DESCRIPTION 


TITL 

NRUNS 

NBLOCK 

NBLOCK 

MSTEP 

CSTRS 

NS UP 


Any alphanumeric description 
of group of runs 

Total number of runs (one run 
corresponds to a unique set 
of input data.) 

Maximum number of blocks to be I 

considered. Crack growth cal- | 

culation ceasos when the number 
of blocks exceeds this munber . 

Block interval for which addi- 
tional data will be printed 
(c.g., 3 would imply that blocks 
3, 6, 9, 12.... etc. would have i 

data printed out). i 

I 

Step interval for which addi- i 

tional data will be printed (in 

the blocks called out above). | 

Constant multiplier for stress I 

inputs. Allows stress spectrum j 

to be varied by changing one 
number only. (e.g., one run with j 
this constant 1 and one run with I 

this constant l«l, will show the | 

effect of varying the stress 1071) . | 

I 

Constant for suppression of re- \ 

tardation in crack growth analysis. 

If zero , retardation is not con- j 

sidcred. If retardation is to be 
considered, constant must be 1. ' 
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CARD 

TYPE 

FIELD 

t«aN|ENCLATURB 

3 

15-18 

NLOAD 

3 

18-22 

NGEOM 

3 

22-26 

NMAT 

3 

27-30 

ITER 

3 

31-40 

PIT 


D^CRIPTION 


Constant to indicato nrhothor now 
load data is to b« input. If saro, 
(or any nunbar not oqual to 1) 
load data will not be read in and 
load data from previous run will be 
used. If it is 1» cord types 4 
and 5 must follow. 

Constant to indicate whether new 
geometry di.ta is to be input. If 
aero (or any number other than 1) 
geometry data will not be read in 
and geometry data from previous run 
will bo used. If it is 1, card 
types 6 and 7 must be read in. 

Constant to indicate whether new 
material data is to be input. If 
zero (or any number other than 1) 
material data will not be read in 
and material data from previous run 
will be used. If it is 1, card 
types 8, 9, 10, and 11, 12, 13, if 
needed must be read in. 

Maximum number of iterations to find 
thickness that produces the desired 
life. May not exceed 10, may be left 
blank. Remaining items on card three 
are left blank, if iter equals zero 
(or blank). 

Parameter to control rate of con- 
vergence on iteration to find thick- 
ness. Usually set to exponential 
power in Paris crack growth equation 
for crack growth data. When other 
crack growth equation is used, app- 
roximate value of a ’’Paris exponent” 
will be sufficient. Nfcist always be 
greater than 1 . 


3 40-50 BLIFE Desired life in blocks. 
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; AHI> 
fYP;: 

FIELD 

:lat«b 

DUj •‘UPTLKI 

4 

1-4 

Niirjp 

Nunbor of at<?ps in load block*. 

4 

4-8 

IR 

Zoi'o If input format includes 
minimum 1 if input 

format inclu Ics stress ratio. 

4 

8-18 

SI>XN 

Limit strrsiv for additional end 
of life rj-ticrmliutlonal, Failare 
due to limit load does not ter- 
mi.iute craiik growth calculation. 

5 

I-IO 

SMAX 

Maximum stress. 

S 

11-20 

SMIN 

Minimum stress if IR 0, stress 
ratio of IK « 1. 

5 

21-30 

UNIT 

Number of cycles or alternate 
rate va'’iaMe. 

5 

31-34 

TYPE 

Mii'*rtvlal property data typo to 
)>e used. 

6 

1-4 

KTYPO 

Initial stress intensity type. 

7 

1-10 

V/ 

Plate width. 

7 

11-20 

TM 

Plate thickness 

7 

21-30 

00 

Initial half surface length for 
oat * t ir )ij 1 • crat'k or center crack. 
( .j 1 , 1 1 t 1 in E i')')-73 

des 1 Lpt ion >ji compact specimen. 

7 

30-40 

• 

AO 

Crack dep^h for a part through 
cra^.k Mav oe left zero for 

through cra<.ks. 

« 

1-4 

NJ 

Number of material property types. 

9 

0 

1 

SIGYS 

Yield stress. 

9 

11-14 

NEQ 

Eountlon to be used for crack 
growth 1 « Col lipriestoEhrct, 

2 * Hans, 'i c Forman. 
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CARD 

TYPE 

PIBLO 

NOMKNCLATURB 

DESCRIPTION 

9 

1S.18 

NRET 

Model to be used for retardation 
0 ■ none, 1 ■ Willonberg, 3 ■ 
Wheeler, 3 ■ Grunuan Closure 
(not debugged). 

9 

18-23 

NDUP 

Constant to indicate whether 
crack growth properties are the 
same in depth and surface dir- 
ections. If constant ■ 1, they 
are not and two sets of D's and 
CR's (see cards 10-13) must be 
input . 

9 

33-33 

KCRC 

Critical stress intensity in 
surface direction (upper cutoff). 

9 

33-42 

KUC 

Threshold stress Intensity in 
surface direction (lower cutoff). 

9 

43-52 

KCRA 

Critical stress intensity in depth 
direction. Need not be input if 
crack is a through crack. 

9 

53-62 

NOA 

Threshold stress intensity in 
depth direction. Need not be 

input if crack is a through crack. 

10 

1-10 

D(1,I,J) 

Constants in crack growth equations 
surface direction. Data is read 
in until a zero is reached. See 
text for description of constants. 

11 

1-10 

• 

I>(2,I,J) 

Constants in crack growth equation 
depth direction. Data is read in 
until a zero is reached. If NDUl* / 
1, card is not used. See text for 
description of constants. 

12 

1-10 

CR(l,I,J) 

Constants in retardation equation - 
surface direction. Data is read in 
until a zero is reached. See text 
for description of constants. 

13 

1-10 

CR(2,I,J) 

Constants in retardation equation - 
depth direction. Data is read in 


until a zero is reached. If NIHJP ^ 
1, card is not uscxl. See text for 
description of constants. 
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CAkD FIELD »VHEN INPUT DATA RBQUIRBO 

TYPE 


1 

1-40 

Once for each session . 



2 

1-16 

Once for each session. 



3 

1-22 

Once for each run. 



3 

22-50 

Whenever an iteration on thickness 
meet design life is desired. 

to 

4 

1-18 

Whenever NLOAD ■ 1 . 



S 

1-34 

NSTEP times vAxcn NLOAD ■ 1. 



6 

1-4 

Whenever NCtiOM ■ 1. 



7 

1-40 

Whoncn/er NGEOM ■ 1. 



8 

1-4 

Whenever NMAT ■ 1 . 



9 

1-62 

NJ times, when NMAT ■ 1. 



10 

1-10 

NJ times, when NMAT * 1. 



11 

I-IO 

NJ times, when NOUP ■ 1 and .NHXT 

■ 

1. 

12 

I-IO 

NJ times, when NMi\T ■ 1 and NRET 

^ 0. 

13 

1-10 

NJ times, when NMAT ■ 1 and NRET 
NDUP ■ 1. 

¥ 

0 and 
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CONTHOL OP OUTPUT 


In addition to information doacribing tho input data and nn 
fracturo, tho computor program output contista of crack longhtOi 
atroaa intonaity factora, crack growth ratoa, and cyclo counts. 
Thoao aro printed out for the first and last cycle in the step. 

This data is always printed for each step of the first block 
encountered in any of the growth nodules (PTCGRW, TRANS or TCCROW) 
Thia data may also bo printed out for additional blocks and steps 
as desired by the user. These additional blocks and steps are con- 
trolled by specifying the increment for blocks and stops for which 
pri. 1 t out will bo made NSiLOCK AND ^STCP respectively. Thus if every 
step in every other block is wiinted, MM/x;K is sot equal to 2 and 
MSTEP is set equal to 1. 
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INPUT FOR ADDITIONAL RUNS 


A full «et of input data io not necessary for each additional run. 

The use of a stress multiplier constant (CSTRS) allows all the stress 
to be varied by a constant percentage without inputting any addi- 
tional input cards other than card type 3. If the stresses are to 
be used directly as they are on card type 5, CSTRS is input as 1. 

The input constant NSUP allows retardation to be suppressed on sub- 
sequent runs. That is if a run is made that considers retardation, 
the following run will perform the same analysis without retardation if 
NSUP ■ 0. When retardation is considered NSUP must be set equal to 1. 
Obviously the order of running the cases must be retarded, followed 
by unretarded. 

In order to control whether loading data, geometry data, or material 
properties are to be read in for a particular run, the constants NLOAD, 
NGEOM, NMAT must be input. If data is to be read in, the appropriate 
constant must be 1, it it is not to be read in, the appropriate con- 
stant is 1 . If data is not read in, data from the previous run is 
used. Obviously, for the first run NIXJAD, NGEOM AND NMAT must all 
be 1. 
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ITERATION OM THICKNESS 


For a given design life (in blocks) the computer program will 
search for the thickness which will meet that life requirement. 

The number of iterations attempted is input by the user. The 
maximum that this may bo is ten. This computer program ceases its 
search when the allawed number of iterations is exceeded or the 
computed life lies between lOOH and 105.1 of the design life. Life 
is arbitrarily defined as the .number of blocks completed plus the 
number of steps completed/total number of steps. 

In addition to inputting the number of iterations and the design life, 
an exponent which will control the rate of convergence to the correct 
thickness must be inserted. iVhen tne Paris equation is used with zero 
threshold for a through the thickness crack, the use of the exponent 
of the Paris equation for the convergence parameter should result in 
a convergence to the correct solution in a single cycle. Any constant 
equal to or greater than the "Paris coefficient" should insure 
convergence. 
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CRACK GROWTH RATE DATA 


Th« following data for use in the Collipr iest/Ehret equation 
(NEQ ■ 1) is typical data for the materials listed. The data 


is included for 

example purposes 

only and 

caution is advised with 

regard to design implications of 

the data 

presented. The 

crack 

growth rate is 

in in/cycle and the stress 

intensity factor 

is in 

Ksi ^in. 

Material 

D(NC,1,J) 

D(NC,2,J) 

D(NC,3,J) 

D(NC,4,J) 


Coefficient 

Relates to Upper 

Lower 

2024-T851 

1.6 X lO"^ 

Midrange 

Slope 

3.45 

Asymptote 

38.0 

Asymptote 

3.4 

2124-T851 

3.3 X 10*^^ 

4.0 

31.0 

3.5 

2219-T87 

2.2 X 10"^ 

3.3 

40.0 

5.5 

(70 F) 

2219-T87 

8.9 X 10"^^ 

4.82 

50.0 

5.5 

(-320 F) 

7075-T6 

4.4 X 10 ® 

2.53 

33.0 

3.0 

7075-T76 

6.3 X 10“^ 

3.0 

30.0 

3.0 

7075-T73 

1.07 X 10’® 

2.67 

40.0 

3.5 

Ti-6A1-4VSTA 

6.8 X 10"^® 

3.3 

50.0 

7.0 

Ti-6AL-4V 

5.7 X 10"'° 

3.18 

84.0 

6.0 

Annealed 

Inconnel 718 

4.0 X 10"^° 

2.7 

115.0 

15.0 

(STA) 

D6AC 

7.5 X 10’^° 

2,74 

90.0 

6.0 
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CRACK GRg^TH ANALYSIS EXAMPLE 


Sample input data for a series of two runs is given on page 29. 

The first run is a part through crack subjected to a series of 
two loads of 1000 cycles each. The material properties associated 
with each load level are different so as to model the effect of 
temperature variations. The second run is the same except that 
the initial crack configuration is a through the thickness center 
cracked panel. No ouput beyond the normally supplied out is re* 
quested for these runs. The output is show'i on pages 30 through 
33. 


A second example co.isisting of a series of three runs is presented 
on pages 3 4 through 45 . Page 34 is the input for these three 
runs. For these runs data on every other block was requested. The 
output was therefore rather extensive. All the Information was not 
needed here and therefore only the first and last page of output of 
those runs discussed are reproduced. 

The first run has two loading steps, calls for a retardation model 
(WiLlenborg) and requests a maximum of three Iterations to find a 
thickness compa*^able with the design life. A rather low (compared 
to the •’Paris exponent") convergence exponent of 2 was used. The 
output for the first iteration is shown on pages 36 through 38 and 
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Th« report on the iteration is shown on page 39 . Due to the 

use of the excessively low value of the co.-ivergence ex{>onent, 
the thickness has not converged to the appropriate value. Note, 
however, that the information is still quite useful and that a 
simple hand plot of the result will show that the correct thick- 
ness is about .52 inches. 

In the second run retardation was suppressed and r.o iteration was 
requested. Note that the life for a thickness of .5 inches goes 
from 137 blocks with retardation to 01 blocks without retardation. 

In the third run a step which similatcd a 480 second hold at a con- 
stant load with a resulting sustained load crack growth was Included. 
The crack growth model was assumed to follow a "Paris" format 
(NSQ B 2) with an exponent of 1 and an appropriate constant was 
assumed. As can be seen by the results, the life was reduced 
further under these assumptions. 
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